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Abstract 

It is likely that the observed large-angular-scale anisotropies in the mi- 
crowave background radiation are induced by the cosmological perturbations 
of quantum-mechanical origin. Such perturbations are now placed in squeezed 
vacuum quantum states and, hence, are characterized by large variances of 
their amplitude. The statistical properties of the anisotropies should reflect 
the underlying statistics of the squeezed vacuum quantum states. The theo- 
retical variances for the temperature angular correlation function are derived 
and described quantitatively. It is shown that they are indeed large. Unfortu- 
nately, these large theoretical statistical uncertainties will make the extraction 
of cosmological information from the measured anisotropies a much more dif- 
ficult problem than we wanted it to be. 
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I. INTRODUCTION 



The line of reasoning in this paper is as follows. 

We see the anisotropies in the microwave background radiation at the largest angular 
scales [1]. Observers convincingly argue that this is a genuine cosmological effect. 

If the large-angular-scalc anisotropy in the microwave background is really produced by 
cosmological perturbations (density perturbations, rotational perturbations, gravitational 
waves), then their today's wavelengths are of the order and longer than today's Hubble 
radius In- Strictly speaking, all wavelengths give contributions to the anisotropy at every 
given angular scale. But if the spectrum of the perturbations is not excessively "red" or 
"blue", the dominant contribution is provided by wavelengths indicated above. In fact, 
the major contribution to the quadrupole anisotropy is provided by wavelengths somewhat 
longer than Ih- 

In the expanding Universe, the wavelengths of perturbations increase in proportion to 
the cosmological scale factor. The wavelengths that are longer than some length scale 
today have always been longer than that scale in the past. Moreover, the wavelengths of the 
perturbations of our interest are much longer than the Hubble radius defined at the previous 
times, when one goes back in time up to the era of primordial nucleosynthesis — the earliest 
era of which we have observational data. It is hard to imagine (although it does not seem 
to be logically impossible) that cosmological perturbations of our interest, with such long 
wavelengths, could have been generated by local physical processes during the interval of 
time between the era of primordial nucleosynthesis and now. We are bound to conclude that 
these perturbations were generated in the very early Universe, before the era of primordial 
nucleosynthesis. There is still 80 orders of magnitude, in terms of energy density, to go from 
the era of primordial nucleosynthesis to the Planck era; a lot of things could have happened 
in between. 

The law of evolution of the very early Universe is not known, but it is likely that it 
could have been significantly different from the law of expansion of the radiation-dominated 
Universe. If so, some amount of cosmological perturbations must have been generated 
quantum-mechanically, as a result of parametric interaction of the quantized perturbations 
with strong variable gravitational field of the very early Universe. Gravitational waves have 
been generated inevitably, while density and rotational perturbations — if we were lucky; 
see [2] and references therein. (If the cosmological scale factor has always been the one 
of the radiation-dominated Universe, we must stop here, because the parametric coupling 
vanishes in this case, and cosmological perturbations cannot be amplified classically and 
cannot be generated quantum-mechanically.) The amount and spectrum of the generated 
perturbations depend on the law of evolution of the very early Universe (the strength and 
variability of the gravitational pump field), and this is how we can learn about what was 
going on there. In particular, the law of evolution of the very early Universe could have 
been of inflationary type. 

If the cosmological perturbations were generated quantum-mechanically, they should 
now be placed in the squeezed vacuum quantum states [3] (for an introduction to squeezed 
states see, for example, Ref. [4,5] and the pioneering works quoted there). Squeezing of 
cosmological perturbations might have degraded by now at short wavelengths but should 
survive at long wavelengths, especially in the case of gravitational waves. 
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Now, the squeezed vacuum quantum states can only be squeezed in the variances of phase 
which unavoidably means the increased variances in amplitude. The statistical properties of 
the squeezed vacuum quantum states are significantly different from the statistical properties 
of the "most classical" quantum states — coherent states. This is well illustrated by the fact 
that the variance of the number of quanta in a strongly squeezed vacuum quantum state is 
much larger than the variance of the number of quanta in the coherent state with the same 
mean number of quanta (N), (N) ^ 1. For a squeezed vacuum state the variance is <iV^> 
- <7V>2= 2 <A^> {<N> +1), while for a coherent state it is <N'^> - <N>'^=<N>. 

The statistical properties of squeezed cosmological perturbations will inevitably be re- 
flected in statistical properties of the microwave background anisotropies caused by them. 
Squeezing is a phase-sensitive phenomenon, and to fully extract its properties the quantum 
optics experimenters use the phase-sensitive detecting techniques based on a local oscillator. 
In cosmology, we are very far from being able to build a local oscillator, except of maybe, 
in the distant future, for short gravitational waves. Besides, in our study of the microwave 
background anisotropies, we are interested in so long-wavelength perturbations that it would 
take billions of years to wait for seeing the time dependent oscillations of variances in the 
quadrature components of the perturbation field. On the other hand, the amount of cosmo- 
logical squeezing is enormously greater than what is achieved in quantum optics laboratory 
experiments. In cosmology, we can only rely on the phase-insensitive, direct detection. One 
can expect that the underlying large variances of the amplitude of cosmological perturba- 
tions should result in large statistical deviations from the mean values for the microwave 
anisotropies. 

A detailed study and proof of this statement is the purpose of this paper. 

At this point it is necessarry to say that the quantum-mechanical generating mechanism 
has become popular in the context of the inflationary hypothesis. However, the inflationary 
literature associates the explanation of the phenomenon with such things as ambiguity in 
the choice of time in the Dc Sitter universe, temperature, tremendous inflation of scales, 
and so on. The basic concepts arc adjusted accordingly. Instead of amplification^ with the 
emphasis on a nonvanishing parametric coupling, increase of amplitude at the expense of 
energy of the pump held, quantum-mechanical generation of waves (particles) in strictly 
correlated pairs, etc., inflationary literature speaks about magnification, with the emphasis 
on "stretching the waves" and "crossing the horizons". Inflationists did not get puzzled 
with their "standard" formula for density perturbations which states that one can produce 
arbitrarily large amount of density perturbations by practically doing nothing. 

The "standard" formula relates the amplitude of density perturbations today with the 
values of the scalar field during infiation. Let us consider, for definiteness, perturbations of 
the matter density Sp/p with today's wavelengths of the order of today's Hubble radius Ih- 
The "standard" formula says that 



Sp 
P 



where the right hand side of this formula is supposed to be evaluated at the time ti when 
the wavelengths of our interest were "crossing the horizon" during inflation. Let us agree 
with the so-called "slow-roll" approximation and assume that the Hubble parameter H was 
almost constant during that epoch, H -C H"^- Let us take the numerical value of H during 
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that epoch at the level, say, 10 orders of magnitude smaller than the Planck value of H. For 
quantum-mechanically generated gravitational waves, it would result in the today's ampli- 
tude h ~ 10~^° and the induced anisotropics of the microwave background 5T/T pa 10~^° 
which are much lower than what is currently discussed in the experiment. However, for den- 
sity perturbations, according to the "standard" formula, the situation is totally different. 
Without changing anything in the curvature of the space-time responsible for the gener- 
ating process (that is, leaving H almost constant and at the same numerical level), but 
simply sending (j)(ti) to zero (which corresponds, due to Einstein equations, to sending H 
to zero, i.e., making the "slow- roll" approximation better and better, making the expansion 
law closer and closer to the De Sitter expansion) one produces arbitrarily large {5p/p)\H- 
Inflationists love to stress that the De Sitter gravitational pump field generates perturba- 
tions with the Harrison-Zeldovich (scale-invariant, fiat) spectrum. What they do not stress 
is that, according to their "standard" formula for density perturbations, the amplitudes of 
the scale-invariant spectrum are infinite, and the amplitudes of the almost scale-invariant 
spectrum are almost infinite. Instead of blaming their own formula, infiationists blame the 
scalar field potentials. This formula is the reason for rejecting certain scalar field potentials 
on the grounds that they generate "too much" of density perturbations, for claims that the 
contribution of gravitational waves to 6T/T is "negligibly small" in the limit of the De Sit- 
ter expansion, and even for claims about copious production of black holes during infiation. 
The "standard" formula has been recently reiterated [6] . It has been formulated, essentially, 
as the following "standard result": "... we see that the scalar perturbations can be very 
strongly amplified" [the increase of numerical value from (almost) zero to (almost) infinity] 
"in the course of the transition" [the instantaneous change of the cosmological scale factor 
from one power-law behavior to another power- law behavior]. The authors of the paper [6] 
assure the trusting reader: "We think that there is nothing strange about this ...". 

The "standard" infiationary formula is unacceptable physically, and it does not come 
as a surprise that it is incorrect mathematically. It may follow from reach imagination, 
but it does not follow from equations. (For technical details see [2,7] especially Exer- 
cises 1, 2, 3 in Ref. [7].) According to the calculations of Ref. [2], the contribution of the 
quantum-mechanically generated gravitational waves to the large-angular-scale anisotropy is 
greater (even in the limit of the De Sitter expansion) than the contribution of the quantum- 
mechanically generated density perturbations. However, in this paper, we will not discuss 
any longer the relative contributions to 6T/T supplied by cosmological perturbations of dif- 
ferent nature (density perturbations, rotational perturbations, or gravitational waves). We 
will concentrate on the consequences of their common origin — quantum mechanics and 
squeezing. Our discussion will be equally well applicable to the perturbations of any nature 
if they have the same origin. 
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II. THE GENERAL EQUATIONS FOR QUANTIZED COSMOLOGICAL 

PERTURBATIONS 



Here we will briefly summarize some basic information about quantized cosmological 
perturbations (see [2,7] and references therein). 

The metric of the homogeneous isotropic universe can be written in the form 



where is the spatial metric. For reasons of simplicity, we will be considering only spatially 
flat universes, that is — 5ij. 

Following Lifshitz, it is convenient to write the perturbed metric in the form 



where hij are functions of 77-time and spatial coordinates. By writing the perturbed metric 
in this form we do not lose anything in the physical content of the problem, but we gain 
considerably in the mathematical tractability of the perturbed Einstein equations. The 
one who is interested in solving equations (not just in discussing them) will certainly be 
interested in a simpler form of equations, even if only as "an insurance against mistakes" 
(the term is borrowed from [6]) when searching for solutions. Those who prefer "gauge- 
invariant formalisms" are welcome to take the found solution and compute with its help 
whichever "gauge-invariant quantity" they like. These quantities, being "gauge- invariant" , 
have the same values in all "gauges". 

The components hij of the perturbed gravitational field can be classified in terms of 
scalar, vector, and tensor eigenfunctions of the Laplace differential operator. The compo- 
nents of the perturbed energy-momentum tensor can also be classified in the same manner. 
After that, the linearized Einstein equations reduce to a set of ordinary differential equa- 
tions, separately for scalar (density perturbations), vector (rotational perturbations), and 
tensor (gravitational waves) parts. 

The number of independent unknown functions of time that can potentially be present 
(on grounds of the classification scheme) in the perturbed Einstein equations is always 
greater than the number of independent equations. It is 6 functions and 4 equations for 
density perturbations, 3 functions and 2 equations for rotational perturbations, and 2 func- 
tions and 1 equation for gravitational waves. In order to make the system of equations closed, 
it is necessary to say something about the perturbed components of the energy-momentum 
tensor or to specify from the very beginning the form of the energy-momentum tensor. The 
popular choices are perfect fluids and scalar fields. Even for gravitational waves, it is not a 
totally trivial question what their definition is (see, for example, Ref. [8]). However, after 
everything is being set, and as soon as the scale factor a{ri) (the background solution) is 
known, the general solution to the perturbed equations can be found. In practice, exact so- 
lutions are being found piecewise, at the intervals of evolution where the energy-momentum 
tensor has simple prescribed forms. 

We can now write the quantum-mechanical operator for the perturbations of the gravi- 
tational field hij in the following universal form: 



ds"^ — —o?{r]){drf' — dx^ dx-') 



(1) 



ds'^ = -a'^{r])[drj'^ - {Sij + hij)dx^ dx^] 



(2) 
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(3) 
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We will start the explanation of Eq. (3) from the polarization tensors Pij. Let us intro- 
duce, in addition to the unit wave- vector n/n, two more unit vectors Zj, m^, orthogonal to 
each other and to n: 

— — (sin 9 cos (p, sin 9 sin 0, cos 9) , k — (sin 0, — cos 0, 0) , 

rrii = ±(cos 9 cos 0, cos 9 sin 0, — sin 9) , (4) 

+ for ^ < f , - for e>l. 

The two independent polarization tensors, s = 1,2, for each class of perturbations, can 
be written as follows. For gravitational waves: 

1 2 

Pij{n) = [klj - mirrij) , (n) = [lirrij + Ijirii) . 

For rotational perturbations: 

1 1 2 1 

Pjj(n) = -{kuj + IjUi) , Pjj(n) = -{miUj + rrijUi) . 

For density perturbations: 

s . s s 

The polarization tensors of each class satisfy the conditions Pijp^^ = 26,^sr, Pij{—n) = Pij{n). 

In practical handling of the density perturbations it proves convenient to use sometimes, in 

1 

addition to the sca/ar polarization component Pij, the longitudinal-longitudinal component 

2 

(proportional to niUj) instead of Pij. The explicit functional dependence of the polarization 
tensors is needed for the calculation of various angular correlation functions. 

The evolution of the creation and annihilation operators Cn(?7), c^{ri), for each class of 
perturbations and for each polarization state, is defined by the Heisenberg equations of 
motion: 

^^-^[cM,H], ^^-^[cl{v),H]. (5) 

The dynamical content of the problem is determined by the Hamiltonian H. Its form depends 
on the class of perturbations and additional assumptions about the energy-momentum tensor 
which we have to make, as was discussed above. 

Under the simplest assumptions about gravitational waves (waves interact only with the 
background gravitational field, there is no anisotropic material sources) the Hamiltonian for 
each polarization component takes on the form 

H = nclcn + ncL„c_„ + 2c7(r^)c^^cl„ + 2a* {rj)cnC-n (6) 

where the couphng function a(r]) is a{r]) = f ^• 

For rotational perturbations, assuming that the primeval matter is capable of supporting 
torque oscillations, assuming that the oscillations are minimally coupled to gravity, and 
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assuming that the torsional velocity of sound is equal to the velocity of light, the Hamiltonian 
for each polarization component reduces to exactly the same form (6) with the same coupling 
function a{r]). 

For density perturbations, we consider specifically a minimally coupled scalar field with 
arbitrary scalar field potential as a model for matter in the very early Universe, and perfect 
fluids at the later eras. The quantization is based on the scalar polarization component (the 
function of time responsible for another polarization state is not independent). There is 
only one independent sort of creation and annihilation operators in this case. The operators 

Cn(^), (^niv) expressible in terms of the operators dnit]), dl^{ri) for which the Hamiltonian 
has again the same form (6) but with the coupling function a{r]) = f-^^^-, where 

For density perturbations, it is the operators dn{rj), dl^{rj) that participate in Eqs. (5), (6). 

Now, let us turn to the constant C in Eq. (3). Its value is determined by the normalization 
of the field of each class to the "half of the quantum in each mode" . Under the assumptions 
listed above, one derives C = yl6n Ipi for gravitational waves, C = y/32n Ipi for rotational 
perturbations, and C — -\/247r Ipi for density perturbations, where Ipi is the Planck length, 

ipi = {Gh/c^y/^ 

The form of the Hamiltonian (6) dictates the form of the solution (Bogoliubov transfor- 
mation) to Eq. (5): 

cU^)=<(^)d(0)+<(^)c-n(0) (7) 

where Cn(0), cJj(O) are the initial values of the operators taken long before the interac- 
tion with the pump field became important {a{r])/n — > 0) and which define the vacuum 
state Cn(0)|0) = 0. The complex functions Un{r)), Vn{rj) obey coupled first-order differential 
equations following from Eq. (4) and satisfy the condition \un\'^ — |f„p = 1 which guaran- 
tees that the commutator relationship [cn(0), cJjj(O)] = 5^(n — m) is satisfied at all times, 
[cn(?7), cjjj(?7)] = (5^(n — m). If one introduces the function jiniv) = '^n{v) + '^niv)^ o^ie recov- 
ers from the equations for Un{r]), Vn{r]) the classical equations of motion. For gravitational 
waves: 

/^n = . (8) 

For rotational perturbations: 

/^n = . (9) 

where Vt is the torsional velocity of sound which we assumed above to be c. For the scalar 
field density perturbations: 

/^n = . (10) 



n — 



a 



2 («V^)" 

n — 
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In the Schrodinger picture, the initial vacuum quantum state |0n) |0_n) evolves into a 
two-mode squeezed vacuum quantum state. In our problem, each of the two-mode squeezed 
vacuum quantum states is a product of two identical one-mode squeezed vacuum quantum 
states which correspond to the decomposition of the real field hij over real spatial harmonics 
sinnx and cosnx. In the Heisenberg picture, the initial vacuum quantum state does not 
evolve in time and is the same now. 

By using Eq. (7) one can present the field (3) in the form 



where the functions hniv) hniv) — '^I'^riiv) + '^niv)]- gravitational waves and 

s s s s 

rotational perturbations, the functions hn are simply hn = t^n/ct' where are solutions to 
Eqs. (8), (9) with appropriate initial conditions. For density perturbations, the functions 

s 

hn are derivable from solutions to Eq. (10) in accord with the relationship between c and d 

1 2 it 2^ 

operators. Besides, for density perturbations, we should regard c^{0) = Cn(0), c^{0) = Cn(0) 
in Eq. (11). In all cases, for a given cosmological model, that is for a model in which the 

s 

scale factor a{ri) is known from the very early times and up to now, the functions hn can be 
found from the classical equations of motion with appropriate initial conditions. 

It follows from Eq. (11) that the mean quantum-mechanical value of the field h^j is 
zero at every spatial point and at every moment of time, (0|/ijj |0) = 0. One can also 
calculate variances of the field, that is the expectation values of its quadratic combinations. 
One useful quantity is hijk'K By manipulating with the product of two expressions (11), 
using the summation properties of the polarization tensors, and remembering that the only 
nonvanishing correlation function is 



one can derive the formula 



|c„(0)c„,(0)|0) = 5,,,53(n-n') 



(0|%(7y,x)/i^^(r?,x)|0) = — / nY.\hn{n)? dn 



(12) 



Equation (12) shows that the variance is independent of the spatial point x but does depend 
on time. 

The expression under the integral in formulas such as Eq. (12) is usually called the power 
spectrum (in this case, it is the power spectrum of the quantity hijh^^): 



Pin)^^nj2\hn{v)\' 



(13) 



s=l 



In cosmology, it is common to use the power spectrum defined in terms of the logarithmic 
frequency interval, that is the function 



PyAn) 



(14) 
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{Z from Zeldovich). We are mostly interested in the power spectrum of cosmological per- 
turbations in the present Universe, at the matter-dominated stage. This spectrum is never 
smooth as a function of frequency (wave-number) n. Squeezing and associated standing 
wave pattern of the field make the spectrum an oscillating function of n for each moment of 
time. However, the spectrum is smooth for sufficiently long waves. At a given moment of 
time, it applies to all perturbations whose wavelengths are of the order and longer than the 
Hubble radius defined at that time. Moreover, the smooth part of the spectrum is power- 
law dependent on n if the scale factor a (77) of the very early Universe (the pump field) was 
power-law dependent on ?7-time. 

Let us assume that the scale factor at the initial stage of expansion was 

a{v) = lo\v\'^^ (15) 

where and /3 are constants. If the evolution is governed by a scalar field, the Einstein 
equations require the constant /? to be /3 < —2. The value f3 = —2 corresponds to the 
De Sitter expansion. At later times, the scale factor changed to the laws of the radiation- 

s 

dominated and matter-dominated universes. From solutions for hiv) traced up to the matter- 
dominated stage, one can find 

E \L{V)\' - r^"""'^' and Pzin) ~ f . 

It is convenient to introduce the characteristic amplitude h{n) of the metric perturbations 
defining this amplitude as the standard deviation (square root of variance) of the perturbed 
gravitational field per logarithmic frequency interval. In the long-wavelength limit under 
discussion, this quantity is universally expressed (both, for gravitational waves and density 
perturbations) by the formula [20] 

Mn)~^n^+' . (16) 

Note that the functional form of h{n) is the same for gravitational waves and density per- 
turbations, the difference is in the numerical coefficient (omitted in this discussion) which 
is somewhat in favor of gravitational waves [2,7]. The numerical level of h{n) is mainly 
controlled by the constant /q- 

The spectra of other quantities can be found in the same manner. For instance, in case 
of density perturbations, one can derive the spectrum of perturbations in the matter density 
Sp/p. Since the relationship between Sp/p and the metric perturbations involves the factor 
(nrj)'^ and, hence, involves two extra powers of n, —{n) ~ rv^h{n), one finds 



p p Jo n 

where P|(n) ~ {lli/ll)n^^^^^^ and 



^(n)~^n^+^ . (17) 

P In 
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It follows from Eq. (16) that h{n) is independent of n if f3 = —2. This independence 
corresponds to the original Zeldovich's definition of the "flat" spectrum: all waves enter 
the Hubble radius with the same amplitude. If the gravitational field perturbations h{n), 
regardless of their wavelength, have the same amplitude upon entering the Hubble radius, 
the matter density perturbations —{n) do also have the same amplitude (the extra factor 
(nr))^ is of the order of 1 at the time when a given wave n enters the Hubble radius). For 
models of the very early Universe governed by a scalar field, the spectral index /3 + 2 in 
Eq. (16) can never be positive. 

Formula (16) and the associated formula (17) should be compared with the "standard" 
inflationary formula which requires that the amplitudes of density perturbations taken at 
the time of entering the Hubble radius should go to inflnity in the hmit of the De Sitter 
inflation, f3 —2. It should also be noted that a "disgusting convention" (the term is 
borrowed from Ref. [9]) is often being used according to which one and the same Harrison- 
Zeldovich spectrum is described by the spectral index = for gravitational waves and 
by the spectral index Ug — 1 for density perturbations. Of course, there is no need in this 
convention. In both cases, the metric perturbations with the Harrison-Zeldovich spectrum 
are described by the same spectral index (zero), see Eq. (16). 

We will flnish this section with a short discussion of coherent states. There is no natural 
mechanism for the generation of cosmological perturbations in coherent states, but if there 
were one it would be reflected in many parts of the theory. The interaction part of the 
Hamiltonian (6) would be linear (not quadratic) in the creation and annihilation operators. 
The analogue of Eq. (7) would read 

c„(r/)=e-^"^cn(0) + a„(r/) 

4(^)=e'""4(0) + <(^) , (18) 

where the complex function q;„(?7) is determined by the coupling function in the Hamiltonian. 
On the position — momentum diagram, the evolution (18) of the field operators corresponds 
to displacing the vacuum state without squeezing while the evolution (7) corresponds to 
squeezing the vacuum state without displacing. In terminology of mechanics, coherent states 
are produced by a force acting on the oscillator while squeezed vacuum states are produced 
by a parametric infiuence. In the coherent states, the mean value of the field is not zero. 
The correlation functions would also be different what would eventually be reflected in the 
differing statistical properties of the field (at least, for some quantities). 
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III. QUANTUM-MECHANICAL EXPECTATION VALUES FOR THE 
MICROWAVE BACKGROUND ANISOTROPIES 



In absence of cosmological perturbations, the temperature of the microwave background 
radiation seen in all directions on the sky would be the same, T. Let us denote a direction 

on the sky by a unit vector e. The presence of cosmological perturbations makes the tem- 
perature seen in the direction e differing from T. The temperature perturbation produced 
by density perturbations or gravitational waves can be described by the formula [10]: 

8T 1 r'^i dh - ■ ■ 

where dhij/dr) is taken along the integration path — e^w, r) — rjn — w, from the event of 
reception w = to the event of emission w = Wi = rjji — tje- The formula for rotational 
perturbations is more complicated than (19) [10], and we will leave rotational perturbations 
aside. 

For quantized cosmological perturbations, the ^(e) becomes a quantum-mechanical op- 
erator. Using Eq. (11) we can write this operator as 

^ C 1 n , r r ^ i ^ 



X 



I s* 



5n(0)/„(«;)e-- + g„(0)/„(«;)e 



(20) 



where 



s 



5 / N 1 rf/lr, 



2n dr] 



Having defined the observable ^(e) and knowing the quantum state |0) we can compute 
various quantum-mechanical expectation values. In the laboratory quantum mechanics, 
the verification of theoretical predictions expressed in terms of the expectation values would 
require experiments on many identical systems. An immediate generalization of this principle 
to cosmology would require speculations about outcomes of experiments performed in "many 
identical universes". Without having access to "many universes" we can only rely on the 
mean (expected) values of the observables and on the probability distribution functions as 
indicators of what is likely or not to be observed in our own single Universe. We will return 
to this point in Sec. IV. 

The expected value of the temperature perturbation to be observed in every fixed direc- 
tion on the sky is zero: 

(0|^(e)|0) = . 

One particular measured temperature map is the result of the measurement performed 
over one particular realization of the random process describing cosmological perturbations 
of quantum-mechanical origin. For this realization, the temperature perturbations may. 
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should, and in fact arc, present. Many measurements will not help (except of reducing the 
instrumental noises) in the sense that they all should give identical results, because the 
timescale of the perturbations under discussion is so enormously larger than an interval of 
time between the experiments. If the COBE's map is correct, we will have to hve with this 
map practically forever. 

Let us now compute the expected angular correlation function for the temperature per- 
turbations seen in two given directions on the sky, ei and 62. This correlation function is 
defined as the mean value for the product of ^(ei) and ^(62): 

i^(ei,e2) = (0|^(e0^(e2)|0) . (21) 
By manipulating with the product of two expressions (20) one can derive the formula 

1 1 rwi rwi foo 

i^(ei,e2) = -C^— — / dw dw ci^^e^"^^!"'-^^'') 
4 (27r)3 Jo io J-00 

X E (P.,(n)elel) (p.,(n)e^e^,) L{^)fni^) ■ (22) 

The next step is the formidable task of taking the integrals over angular variables in 3- 
dimensional wave- vector n space. However, it can be done (see Ref. [11] for gravitational 
waves, Ref. [12] for rotational perturbations, and Ref. [2] for density perturbations). The 
final expression reduces, without making any additional assumptions whatsoever, to the 
form 

00 

K{ei,e2)=K{6)=l% J2 KiPii^^^^) ■ (23) 

We see that the correlation function depends only on the angle 6 between the directions ei, 62 
not directions themselves. The coefficient Ipi is taken from C^, other numerical coefficients 

s 

are included in Ki. The quantities Ki involve the integration of /„(w) over the parameter w 
and the remaining integration over the wave-numbers n. The numerical values of Ki depend 
on a chosen sort of cosmological perturbations and a chosen cosmological model; so far, the 
formula (23) is totally general. Pi{cos6) are the Legendrc polynomials. The lowest multipole 
^min follows automatically from the theory and it turns out to be, not surprisingly, /min = 
for density perturbations, l^am — 2 for gravitational waves (and Zmin = 1 for rotational 
perturbations). For the separation angle 5 — 0, Eq. (23) reduces to the variance of y"(^)) 
that is 

Sjrji crri OO 

(0|^(e) — (e)|0) = i^(0) = /|,, E • (24) 

•' — '•min 

Formula (23) gives the expected value of the observable ^(61)^(62). If the experimenter 
measured this observable in "many universes" and averaged the measured numbers, he/she 
would get the result (23). Moreover, formula (23) says that if the experimenter made the 
measurements at any other pair of directions, but with the same separation angle S, he/she 
would again get, after the averaging over "many universes", the result (23). Without having 
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access to "many universes", we can ask what is the theoretical standard deviation of the 
quantity (61)^(62). (In practice, for deriving K{6) we need a kind of ergodic hypothesis 
allowing to replace the averaging over "universes" by the averaging over pixels separated by 
a given angle on a single map.) The variance V{ei, 62) of this quantity is, by definition, 

(25) 



V(e„e,) = (0|^(e0^(e2)^(e0^(e2)|0) 



(0|^(ei)^(e2)|0) 



The standard deviation is the square root of this number. 

The calculation of ^(ei, 62) requires us to deal with the product of four expressions (20). 
However, the mean values of the products of four creation and annihilation operators are 
easy to handle. One can show that ^(ei, 62) depends only on the separation angle 5 and 



V{e,,e2)=V{S) 



, 5T , , 6T , , , , 
I ^(ei) — (62)10) 



+ 



, , 5T , , 5T , , , , 
(0| — (e) — (e)|0) 



that is 



ViS) = K\S) + K\0) 



(26) 



(27) 



(28) 



The standard deviation for the observable ^(61)^(62) is 

In a similar fashion one can derive the higher order correlation functions for two directions 
ei, e2 and the correlation functions for larger number of directions, but we will not need 
this information. 

In the limit 6 = Eqs. (25), (26) say that 

2 



(0| 



ST 



|0) 



6T 



1 2 



|0) 



(29) 



The familiar factor 3 relating the fourth-order moment with the square of the second-order 
moment (given that the first-order moment is equal to zero) is the reflection of the underlying 
Gaussian nature of the squeezed vacuum wavefunctions associated with the Hamiltonian (6). 

By examining Eq. (28) one can conclude that for each separation angle 6 the standard 
deviation of the angular correlation function is very big. Even at those separation angles 
at which K(S) vanishes, the standard deviation is as big as the variance for ^(e) itself. 
However, the value of the standard deviation for a given variable is not very informative 
per se, as long as the probability density function for this variable is not known. If the 
probability density function (p.d.f.) were normal, we could say that the probability to find 
a result outside of la interval is 32%. Without knowing the p.d.f. we could resort to the 
Chebyshev inequality, but it would only tell us that this probability is less than 1. To 
get more information about possible deviation of the angular correlation function from its 
mean values we will consider in Sec. IV a classical random model which will reproduce the 
expectation values calculated above and will allow us to construct the p.d.f. for the variable 
^(ei) ^(62). On the other hand, the quantum-mechanical calculations of this Section will 
shed light on the classical model. As is known, "quantum mechanics helps us understand 
classical mechanics", see on this subject a paper of Zeldovich signed by the pseudonym 
Paradoksov [13]. 
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IV. CLASSICAL MODEL FOR THE STATISTICS OF THE MICROWAVE 

BACKGROUND ANISOTROPIES 



A distribution of the microwave background temperature over the sky is a real function of 
the angular coordinates. Assuming that 5T/T is a sufficiently smooth function on a sphere, 
one can expand it over the set of orthonormal complex spherical harmonics Yi^{6, 0) [14]: 

5T °° ' 

I7r(e) = E E hmYUe) + a^^Y^^ie)] . (30) 

1=0 m=-l 

We want to formulate a statistical hypothesis about the coefficients aim, so it is better to 
write them first in terms of real (r) and imaginary (i) components: 

^im — Yfm + iYil^ , Y*^ — Yf^ — iYil^ , 



oo I 

T 



'^^ e) = 2 E E i^LYM - aLYM] ■ (31) 

/=0 m=-l 

Our statistical hypothesis is as follows: (i) all members of the set of random variables 
{a[^, al^} are statistically independent, (ii) each individual variable is normally distributed 
and has a zero mean, (iii) all variables with the same index I have the same standard 
deviation ai. All said is expressed by the probability density function (p.d.f.) for individual 
variables: 



1 ("Fm)" 1 (°L)^ 

/KJ = ^^e~^, /(aL) = ^^e"^, (32) 

and by the p.d.f. for the entire set of variables, which is simply a product of all p.d.f. 's for 
all individual variables: 

/({«L,«U)=n;,m/(a[J/(aL) ■ (33) 

Having postulated the p.d.f.'s, we can now compute the expectation values of certain 
functions of the random variables. Below, the angular brackets will denote the expectation 
values calculated with the help of the p.d.f. (33), unless other definition is stated. 

Obviously, all linear functions have a zero mean: 

<a\^> =0, <aj^> =0 . (34) 
For quadratic combinations we have 

<«Ui«U2> =0> <«UiaU2> =0. (35) 

All triple products have zero means. Among quartic combinations, only those can survive 
which have four indices (r), or four indices (i), or two indices (r) and two indices {%). Two 
representative expressions are: 
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+ '^li'^l2^hl3^mim3^l2li^m2mA'^ '^Ii'^l2^lil4^mim4^l2l3^m2m3 , (36) 
< ^hmi ^hm2 ^hm3 ^hmA ^ ~ '^h '^'I'a ^mim2 ^hU ^mzmi ■ (37) 

Other quartic combinations can be obtained by the replacement (r) ^ (i) in Eqs. (36), (37) 
(or by permutation of pairs {Im) in case of Eq. (37)). The higher-order correlations can be 
derived in a similar way, but we will not need them. 

In our further calculations related to the random variables ^(e) and ^(ei)^(e2) it is 
easier to deal with the complex coefficients aim-, so we will first translate the above relation- 
ships to them. By using the available information one can derive 

<a/m> =0, <ahmi(l*2m2> ^'^^l^h l2^raim2 ) 
'^Iimi0'l2m2(^l3m3^l4m4 ~ 4:af^af^{Si^igSfnim3^l2l4,^m2m4 + ^Iil4^mim4^l2l3^m3m3) ■ (38) 

The mean values of the complex conjugated quantities are given by the same formulas (38). 
Other nonvanishing quartic combinations can be obtained from the one in Eq. (38) by the 
permutation of pairs (/m). 

Now, even before deriving the p.d.f.'s for the random variables ^(e) and ^(61)^(62), 
we can find some expectation values. It is clear from the definition (30) and Eq. (38) that 

= . (39) 

When calculating the angular correlation function one should remember that 

E YUei)YrJe,) = ^P^cos^) (40) 

m=—l 

(note the origin of the factor 21 + 1 which will accompany us often). By taking the product of 
two expressions (30) and using Eqs. (38), (40) one can find the angular correlation function 

(^(ei)^(e2)^ = -E'^'(2^ + l)^Kcos(5) . (41) 

If the separation angle S is zero, we obtain 

-| 00 

-Y^aU2l + l) . (42) 
^1=0 

[One may notice an incidental fact that the mean value of the random variable defined as 
= Y}m=-i '^imO'lm < >= 2(2/ + l)crf , that is the same expression which enters Eq. (41). 
This may suggest an interpretation of the quantity <af> — <a^>'^ as the variance of the 
multipole moments. It would be an error, see Appendix.] 

Wc can also find the 4th order expectation values. The product of 4 expressions (30) in 
conjunction with Eq. (38) gives 
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'5T ST ST ST \ /ST ST ' 
^(61)^(62)^(61) — (62) ) - ( — (ei) — (62) 

/ST^ ST^ X IST ^ ST^ y 
= (— (ei) — (62)) +(— (e) — (e)\ . 



(43) 



If 5 = 0, it follows from Eq. (43) that 

1 4 



ST 
~T^ 



e) 



ST 



e) 



2\ 2 



(44) 



Up to difference in the meaning of the angular brackets, the formulas (39), (43), (44) 
reproduce the analogous results of the previous Section. Moreover, from comparison of 
Eqs. (23), (24) with Eqs. (41), (42) wc can relate the quantities Ki, derivable from a given 
cosmological model plus perturbations, with the abstract quantities ai. 

We can now engage in our major enterprise — the construction of the p.d.f. for the 
random variable v = ^(61)^(62). We will start from the p.d.f. for the random variable 
z = ^(e). When it is necessary to distinguish directions ei and 62, we will use the notations 
zi and Z2. 

The variable z is a function of the variables {a[„, a;„} whose p.d.f.'s are known, Eqs. (31), 
(32). There exist regular methods (see, for example, an excellent book [15]) allowing to derive 
rigorously the p.d.f. of a function. However, in our case that the function is linear and all 
p.d.f 's are normal, we can partially rely on a guesswork. Combining formulas and guessing 
we can write 



27r (T2 



where 



1 00 
^1=0 



(45) 



(46) 



The p.d.f. (45) certainly leads to Eqs. (39), (42), (44). Moreover, it allows us to say that 
the probability to find z outside of la^ interval is approximately 32%: 

P{\z\ > a^) « 0.32 . 

We now introduce two variables, Zi and Z2, and ask about the p.d.f. in the 2-dimensional 
space (^1, ^2)- Again, partially relying on a guesswork, we find that 



f{Zl,Z2) 



2..V1-P- '"' {~ 2.2(i'-,2) [^i + -2 - 2P.1.2]} (47) 



where 



pal^-^(2l + l)afPi{cosS), \p\ < 1 . 
^1=0 



(48) 
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(Sec Eq. (5.11.1) in Ref. [15]). First, we can check that the marginal distributions are correct. 
For f{zi), one obtains 



r-oo I 
f{Zl) = / f{Zl,Z2)dZ2 = -7= 

j-oo yZir 



4 

2^ 



and one obtains a similar expression for f{z2)- Second, one can check that 

/ 2\ 2 / 2\ 2 / \ 2 

where the angular brackets mean the integration with the p.d.f. (47). These equahties are 
Eqs. (42), (41) which we must have obtained. 

Finally, we shall derive the p.d.f. for the variable v — Z\Z2- We will do this in some detail 
following the prescriptions of [15]. 

Let us introduce the two new variables {z\^v) instead of (-21,-22) according to the trans- 
formation 

V 

Zl = Zi , Z2 = — 
Zl 

The Jacobian of this transformation is J — 1/zi. The p.d.f. f{v) is the result of the 
following integration: 



zl + ^ - 2pv 



V 



2 



z 



dzi . 



The integral over zi can be taken with the help of 3.471.9 from [16]. The resulting p.d.f. 
can be written in the form: 

f{v) = { --'V^ Y'^' ' ^[ (49) 



-^=e^I(^is:o f^T^) , for ^; < 



where Kq is the modified Bessel function of its argument [17]. 

The function f{v) is quite complicated and the distribution is obviously not normal. The 
function f{v) goes to zero for v ±00 and diverges logarithmically at the point v = 0. 
Even a verification of the normalization condition 

/oo 
f{v)dv = 1 (50) 
-00 

is not trivial. However, with the help of 6.621.3, 9.131.1, 9.121.7, 1.624.9 and 1.623.2 
from [16] one can prove the validity of Eq. (50). 

The mean value and the standard deviation of the variable v are known, see Eqs. (41), 
(43): 

1/2 



{v) = pal , a, = [{v') - {v)'] = al4p^l . (51) 

We already knew that the standard deviation is big. We now see again that = at the 
separation angles at which the angular correlation function vanishes, p — 0, and (7„ = V^cr'^ 
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at zero separation angle, p — 1. For other separation angles, a^, lies between these two 
numbers. 

Now that we know the p.d.f., we can assign probabilities to the different ranges of the 
variable v. For instance, we can calculate the probability that the measured v will be found, 
say, outside of the A(T„ interval surrounding the mean value of v, where A is an arbitrary 
fixed number. The probability of our interest is 



To get a quahtativc estimate of the associated theoretical uncertainties for the observable 
V, we will ask a slightly different question. What should the number A be in order to have 
the 0.32 chance of finding v outside the Acr^ interval and, hence, the 0.68 chance to find it 
inside the interval? 

To evaluate the size of the disaster, we will start from the case p = 0. In this case, the 
p.d.f. (49) is symmetric with respect to the origin v — (this is why {v) is zero in this case) 
and 



We want this number to be approximately equal to 0.32. Judging from the Fig. 9.7 in 
Ref. [17], a half of the area under the Kq{x) function is accumulated when integrating from 
approximately a; = 1/2 and up to infinity. This means that A should approximately be equal 
to 1/2 . 

If p 7^ the evaluation of P is more complicated. For p 7^ 0, the function (49) is not 
symmetric with respect to the origin v = 0. It has larger values at positive v^s if p > (this 
is why (v) > in this case) and it has larger values at negative f 's if p < (this is why 
(v) < in this case). The graph of the function e^Ko{x) plotted on Fig. 9.8 in Ref. [17] 
is helpful. A qualitative analysis shows again that A is approximately equal to 1/2. (More 
accurate estimates can of course be reached by numerical methods.) 

At any rate, the ^a^, interval gives approximately the same probability estimates as if 
the distribution (49) were normal. 




(52) 




(53) 
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V. CONCLUSIONS 



A particular cosmological model plus perturbations gives unambiguous predictions with 
regard to the expectation values of the measurable quantities. Differing models give different 
predictions. We want to distinguish them observationally and to learn about physics of the 
very early Universe. However, the quantum-mechanical origin of the cosmological pertiir- 
bations is reflected in the theoretical statistical uncertainties surrounding the expectation 
values. One important measurable quantity is the angular correlation function of the mi- 
crowave background anisotropies. Its mean value at the zero separation angle was denoted 

in this paper. It was shown that the standard deviation for the correlation function is 
very big. The 68% confidence level corresponds, approximately, to |(T^ at the separation 

angles where the correlation function vanishes, to -^cr^ at the zero separation angle, and to 
intermediate numbers for other separation angles. 

The angular correlation function has actually been measured. It is presented at Fig. 3 
in the paper [18]. The authors surround the measured points by a narrow shaded region 
which they address as follows: "The shaded region is the 68% confidence region .... including 
cosmic variance and instrument noise" . It is not quite clear what the authors of Ref. [18] 
(see also Ref. [19]) mean by "cosmic variance", but if they mean the theoretical statistical 
uncertainties, these uncertainties are significantly larger than what is plotted. According 
to the calculations presented above, the half-width of the shaded region should be approxi- 
mately 600 {ij-KY near the points where the correlation function vanishes and approximately 
840 {ixKY near the point marking the zero separation angle. 

The conclusion is a bit disappointing. Apparently, God is telling us something important 
about the very early Universe by exhibiting the microwave background anisotropies, but the 
channel of information is so noisy that it will be hard to understand the message. 
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APPENDIX 



The set of random variables {a[^, a^^} defined by Eqs. (32), (33) lives its own independent 
life regardless of whether or not the variables are considered random coefficients in the 
expansion of some function over spherical harmonics. Being such, it allows introduction of 
new functions and calculation of their expectation values. One interesting variable is defined 
by the equation 

I I I 

m=—l m=—l m=—l 

By using Eq. (38) one can calculate the expectation value of af: 

(af) ^ {21 + l)2af . (55) 

The factor (2/+ 1) reflects the number of independent "degrees of freedom" associated with 
the index /. One can also introduce the variable af and calculate its expectation value: 

(af) = (2/ + l)(/ + l)8af = (af)^^[±|i . (56) 

The difference (af) — {afY is, by definition, the variance of the variable af . Prom Eqs. (56), 
(55) one finds 

{<) - {alf - ^^{air . (57) 
Now, we can notice that the angular correlation function (41) can be written in the form 



oo 



^f^(ei)f^(e2)^ = ^E(«f>^Kcos5) . (58) 



On this ground, there may be a temptation (which is, in fact, followed by many including 

"t"(^i)"t" 



observers) to write the random variable ^(ei) ^-(62) in the form 



00 



f^(ei)f^(e2) = ^E«'^Kcos<5) (59) 

and to interpret Eq. (57) as the variance for the multipole moments of the correlation 
function. One should strongly resist to this temptation. 

Let us show that the definition (59) is incorrect despite the fact that it gives correct 
expectation value (58). It follows from the definition (59) that 

5T ST ST ST 
^(ei)^(e2)^(ei)— (62) 

= T3E«f[^Kcos5)f + — E afPi{cosS)4P,{cosS). (60) 

1=0 l,V={),l^l' 
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Using (56), (58) and remembering that af and af, are statistically independent for I ^ I', 
one can find the expectation value of the quantity (60): 



'6T 6T 5T 6T \ 
jr(ei)^(e2)^(ei) — (62) ) 

= — ^(af)[P,(cos5)]^ + — {a^){4)Pi{cos6)P,{cos6) 



1=0 



-1 00 

^(af)[P,(cos5)]2 + 



^" l,l'=0,l^l 

1 



47r2 



1=0 



47r2 



^(af)P,(cos(5) 
7=0 



-1 00 



47r2 



1=0 



/-(ei)^(e2)) +^,i::^{amPiicosS)r. 



\ T 



47r2^2Z + l 



(61) 



It follows from (61) that the variance of the variable (^1)^(^2) would read (if (59) were 
correct): 



-1 00 



TT 



(62) 



1=0 



This expression should be compared with the correct variance following from Eq. (43): 



1 



1 2 



u=o 



1 

+ — 



E(2^ + 



./=o 



(63) 



Formulas (62), (63) disagree even for 5 — 0, and even in their first, I — term. This shows 
that the ad hoc definition (59) is incorrect. The correct definition of the random variable 
^(ei)^(e2) is the one following from the definition (30) and which we have used in this 
paper. 
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